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We investigate how the microstructure of a colloidal polycrystal influences its linear visco¬ 
elasticity. We use thermosensitive copolymer micelles that arrange in water in a cubic crystalline 
lattice, yielding a colloidal polycrystal. The polycrystal is doped with a small amount of nanopar¬ 
ticles, of size comparable to that of the micelles, which behave as impurities and thus partially 
segregate in the grain boundaries. We show that the shear elastic modulus only depends on the 
packing of the micelles and does not vary neither with the presence of nanoparticles nor with the 
crystal microstructure. By contrast, we find that the loss modulus is strongly affected by the pres¬ 
ence of nanoparticles. A comparison between rheology data and small-angle neutron scattering 
data suggests that the loss modulus is dictated by the total amount of nanoparticles in the grain 
boundaries, which in turn depends on the sample microstructure. 


INTRODUCTION 


Most solid materials in everyday life, as e.g. metals 
or ceramics, are crystals, whose constituents (atoms or 
molecules) are regularly positioned in space. Over the 
past century scientists have reached a good understand¬ 
ing of the electric, thermal, and mechanical properties of 
a defect-free crystal. However, real-life crystalline ma¬ 
terials have defects, which strongly affect their proper¬ 
ties [l|. In particular, crystalline solids such as metals 
and ceramics are usually aggregates of crystalline grains 
that form the micro structure of the material, which plays 
a key role in almost all industrial processes and products. 

Very generally, the mechanical properties of crystalline 
materials strongly depend on defects in the crystalline 
structure. The interplay between defective structure and 
mechanical properties has been extensively investigated 
in metals, and more recently in nano crystalline materials, 
i.e. materials where the size of the grains is smaller than 
typically 100 nm, because of their outstanding chemical 
and mechanical properties, including superplasticity [l- 
[^. Grain boundaries (GBs), the two-dimensional lattice 
defects that separate the different grains, control the bulk 
mechanical properties of polycrystalline materials ii, 
SliMl- In particular, the sliding and the migration of 
GBs play important roles in the plastic deformation and 
the fracture at high temperature. 

How GBs and defects in general impact the mechani¬ 
cal properties of soft materials is far less documented. A 
few exceptions include liquid crystalline materials where 
inclusion of particles has been shown to induce defects 
that strongly modify the linear mechanical properties of 
the composite due to the presence of a stabilized network 
of defects The link between microstructure and 

mechanical properties in colloidal materials was invoked 
in a few previous studies with dispersions of block copoly- 
15l4l8j but was not investigated in detail and we 


are not aware of any experimental work that investigates 
how the visco-elasticity of a colloidal polycrystal depends 
on its microstructure, i.e. on the average size of the crys¬ 


talline grains. A prerequisite for such study is to be able 
to tune the average size of the crystalline grains, which 
in atomic and molecular crystals is typically achieved by 
varying the cooling rate. For colloidal crystals, this is 
not an easy task Pill because the crystallization rate 
is controlled by the volume fraction of colloids, a param¬ 
eter that usually cannot be varied during an experiment, 
unlike temperature T for hard condensed matter systems. 


This difficulty is overcome in a novel colloidal system 
that we have recently developed and that displays analo¬ 
gies with atomic and molecular crystals. First, thanks to 
the thermosensitivity of the colloids, the crystallization 
can be induced at a controlled rate by changing tem¬ 
perature [ii-il- Second, the system comprises small 
amounts of nanoparticles of physico-chemical nature dif¬ 
ferent from that of the colloids forming the crystal, which 


act as impurities in atomic and molecular systems [25 


and as such partially segregate in the grain boundaries 
during the crystallization process [23|, [271, l28| . The mi¬ 
crostructure can be controlled both by the rate at which 
the temperature is varied, and by the amount of nanopar¬ 
ticles. We leverage on these specificities to investigate by 
rheology how the microstructure influences the sample 
visco-elasticity. Gomplementary neutron scattering data 
are used to understand the microscopic origin of the ob¬ 
served mechanical behavior. 


The paper is organized as follows. We first describe the 
materials and experimental methods. We then present 
how the linear visco-elasticity varies with the amount of 
nanoparticles and the temperature rate used to prepare 
the samples. Gomplementary experiments by scattering 
techniques are shown in order to relate our findings on the 
visco-elasticity to the sample microstructure. We finally 
discuss our results in light of the knowledge on atomic 
and molecular crystals, and on other soft materials with 
a cellular structure. 
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MATERIALS AND METHODS 
Experimental samples 

The colloidal crystals have been described else¬ 
where [IMil- They comprise a commercial block- 
copolymer, water and nanoparticles. We use a ther¬ 
mosensitive triblock copolymer of Pluronics type (F108), 
made of a central polypropylene oxide (PPO) block of 
52 monomers flanked at each extremity by two polyethy¬ 
lene oxide (PEO) blocks of 132 monomers each. The 
concentration of F108 in water is set at 34% w/w. The 
affinity of the PPO block for water varies with temper¬ 
ature. At low temperature (T < 6°C), the copolymer 
is fully water-soluble and the sample is a viscous fluid 
(zero-shear viscosity 0.32 Pa.s at T = 4°C). Conversely, 
at room temperature the affinity of PPO for water is 
poor, and FI08 self-assembles into micelles of diameter 
2a = 22 nm that arrange on a face-centered cubic lat¬ 
tice (lattice parameter 32 nm) A small amount 

of nanoparticles (NPs) is added to the F108/water mix¬ 
ture. We used Bindzil silica NPs (kind gift from Eka 
Chemical, sample type 40/130), with an average diam¬ 
eter of 30 nm and a relative polydispersity of 19%, as 
determined by transmission electron microcopy. A few 
complementary experiments are performed with carboxy- 
lated polystyrene nanoparticles of diameter 36 nm pur¬ 
chased from Invitrogen. The volume fraction of NPs, 0, 
is varied in the range (0.5 — 2)%. The sample is pre¬ 
pared in a fluid phase at T ^ 3°C. Crystallization is 
then induced by rising T at a fixed heating rate, with 
T in the range (0.001 — 2)°C/min. We have previously 
shown that the crystallization temperature continuo usly 
increases with the heating rate from 15.5 to 18.5°C [30|. 
Both rheology and scattering experiments are performed 
at T = 23°C, where all samples are fully crystallized. 

For neutron scattering measurements, water {H 2 O) is 
replaced by an equal volume of dent crated water {D 2 O) 
or a mixture of D 2 O and H 2 O, with various ratios, in 
order to match the scattering of the solvent with that of 
the copolymer (resp., with that of the NPs), so that only 
the NPs (resp., only the micelles) should contribute to 
the scattering signal. Samples without NPs are prepared 
with pure D 2 O in order to maximize the contrast between 
the copolymer and the solvent. Samples with NPs are 
prepared using 39/61 H 2 O/D 2 O w/w to contrast-match 
the NPs or 85/15 H 2 O/D 2 O w/w to contrast-match the 
copolymer, as detailed previously [ 2 ^. We have checked 
that changing the solvents does not modify the sample 
structure |29| . 

The NPs act as impurities in the colloidal crystal. As 
such, they are found to partially segregate in the grain 
boundaries upon crystallization. As an illustration, we 
show in fig.[T]a light microscopy image of a sample com¬ 
prising nanoparticles. Thanks to the accumulation of 



Figure 1: Differential interference contrast microscopy im¬ 
age of a sample comprising nanoparticles with a volume 
fraction 0 = 2%, and prepared with a heating ramp T = 
0.001 °C/min. Scale: 20/im. 

NPs in the GBs, the contrast between the crystalline 
grains and the GBs is enhanced and allows for the imag¬ 
ing of the network of GBs that delineate grains with dif¬ 
ferent crystalline orientations. For the sample shown in 
fig.[T](T = 0.001°G/min, 0 = 2 %), the average size of the 
crystalline grains is of the order of 10 /im, but the sample 
microstructure can be tuned by varying the amount 0 of 
NPs, or the heating rate T [23|. The average size of the 
crystallites typically ranges between 8 and 45 jam. 

Using fluorescent polystyrene nanoparticles of size 
comparable to that of the silica NPs used here, and quan¬ 
titative confocal imaging, we have previously shown [ 2 ^ 
that not all NPs are located in the GBs. A partitioning 
coefficient, 7^, can be introduced, defined as the ratio be¬ 
tween the total amount of NPs in the GBs and the total 
amount of NPs inside the grains. We have found [23| that 
V is not constant but rather decreases as the crystallite 
size becomes smaller. Interestingly, V only depends on 
the microstructure and not on what combination of 0 
and T is used to achieve a given crystallite size. 


Experimental techniques 

The linear visco-elasticity of the samples is investigated 
using a stress-controlled rheometer (USD 200, Physica) 
operating in a strain-controlled mode thanks to a feed¬ 
back loop. The sample, initially in a liquid phase at 
T = 3°G, is loaded in a Gouette cell whose temperature 
is also set at 3°G, and is covered with a thin layer of sili¬ 
con oil to prevent evaporation. The sample temperature 
is increased at a fixed temperature rate T using a water 
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circulating thermostat bath. Visco-elasticity measure¬ 
ments are performed by imposing an oscillatory strain 
of amplitude 7 and frequency / and by measuring the 
complex shear modulus, G* = G' + iG" (ssj. 

Additional experiments are performed using small- 
angle neutron scattering. Details of the experimental 
conditions are described elsewhere [ 2 ^. In brief, experi¬ 
ments are run on the beamline PACE at the Laboratoire 
Leon Brillouin, Saclay, France. To cover the range of 
scattering vectors q G [(0.003 — 0.4)A“^], three experi¬ 
mental configurations are used, D = 4.5 m, A = 12 A ; 
D = 4.5 m, A = OA, D = 1 m, A = OA, where D is 
the sample-detector distance and A is the neutron wave¬ 
length. Here, q = 47 rA“^ sin( 6 >/ 2 ) is the magnitude of the 
scattering vector and 6>/2 is the scattering angle. Empty 
cell background subtraction, calibration by light water 
in a 1 mm-thick quartz cell, and absolute determination 
of scattering intensities I{q) per sample volume are per¬ 
formed using standard procedures 29|, . 


RESULTS AND DISCUSSION 

Visco-elasticity of samples without nanoparticles 

We first investigate the linear visco-elasticity of sam¬ 
ples without nanoparticles. Strain sweeps measured at a 
frequency of / = 0.5 Hz and frequency sweeps measured 
in the linear regime (7 = 0 . 1 %) are shown in fig. [ 2 ] for 
samples prepared with different heating rates T. Over¬ 
all, the visco-elasticity is only very weakly dependent 
on T. The 7 dependence of the complex modulus ex¬ 
hibits the hallmarks of the visco-elasticity of soft solids, 
as found for many materials ranging from colloidal gels 
to concentrated emulsions and dry foams [35|. In the 
linear regime, below 7 0 . 2 %, the storage modulus G' 

and the loss modulus G” are independent of the applied 
strain. Above a critical strain of order 10 %, G' oc 7 “^. 
The loss modulus, by contrast, exhibits a well defined 
peak, whose maximum occurs roughly at the cross-over 
of G' and G”, before decreasing proportionally to 7 “^. 
The frequency-dependent complex modulus in the linear 
regime is dominated by the storage modulus, which, in 
the range of frequency investigated [(0.005 — 20) Hz], is 
between one and two orders of magnitude larger than the 
loss modulus, indicating a highly elastic soft material. 
We find that the storage modulus is nearly frequency- 
independent over the whole range of probed frequencies. 
Moreover, the data for G' acquired at different heat¬ 
ing rates T G (0.007 — 2)°C/min perfectly superimpose, 
yielding a frequency- and heating rate-independent elas¬ 
tic plateau. Go = (11610 ±640) Pa (fig. [3]A). We find on 
the other hand that the loss modulus continuously de¬ 
creases from about 1000 Pa down to about 100 Pa as the 
frequency increases from 0.005 to 20 Hz. At relatively 
high frequency (/ > 0.5 Hz) the loss modulus slightly 




/(Hz) 


Figure 2: (Color online). Strain sweeps (A) and frequency 
sweeps (B) for samples without nanoparticles, prepared using 
various different temperature ramps as indicated in the leg¬ 
end. The frequency is 0.5 Hz in the strain sweep experiments 
and the strain amplitude is 0.1% in the frequency sweep ex¬ 
periments. 


varies with the heating ramp. At / = 2 Hz, G” in¬ 
creases from 110 to 215 Pa and finally reaches a plateau 
for T > 0 . 1 °C/min (circles in fig. [3B)- 


Visco-elasticity of samples doped with nanoparticles 

The elastic modulus of samples doped with nanopar¬ 
ticles behaves closely to that of samples without NPs, in 
that G' is both frequency-independent (data not shown) 
and nearly independent of T (fig.[3]A). Moreover, the nu¬ 
merical value of the elastic modulus is insensitive to the 
content of NPs: in the explored range 0 < 0 < 2%, all Go 
vs T data collapse almost perfectly on the same master 
curve as for samples without NPs, as seen in fig. [3K- Fi¬ 
nally, the way the storage modulus depends on 7 in the 
non-linear regime is similar to that for a sample with¬ 
out NPs. This is illustrated by the superimposition of 
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the data with and without NPs for 7 > 3 %, i.e. in the 
regime where G” > G' (fig. [4]A.). As seen in fig. |4l4, this 
also holds for the loss modulus G". Hence, the addition of 
NPs does not alter the elastic properties of our polycrys¬ 
tals in both the linear and non-linear regimes, nor does 
it modify the viscous properties in the large-amplitude, 
non-linear regime. 

By contrast, in the linear regime the loss modulus ex¬ 
hibits significant variations between samples with and 
without NPs. We first discuss the differences between a 
pure polycrystal and samples doped with a fixed amount 
of NPs and prepared at various ramp rates T. Note that 
we describe the data obtained with (j) = 1 %, but the 
main conclusions are robust with respect to the amount 
of added NPs. First, in sharp contrast to the behavior 
of samples without NPs, the strain sweeps performed at 
/ = 0.5 Hz show that, at low 7 , G” strongly increases as 
the heating rate decreases (fig. IU\). We plot in fig. [SB 
G"(/ = 2 Hz) as a function of the temperature ramp 
rate. The loss modulus decreases when the heating rate 
increases and reaches a plateau for T > 0.1°C. This trend 
is robust with respect to the choice of the frequency at 
which the loss modulus is measured: as it can be inferred 
from fig. IIB, a decrease of G" followed by a plateau at 
T > 0.1°C is systematically seen in the whole range of 
accessible frequencies. (An example for data acquired at 
/ = 0.05 Hz will be shown in fig. [ 6 ^ below). Within this 
general trend, only the amplitude of the decrease of G" 
with T changes, the largest drop of the loss modulus be¬ 
ing observed at low frequency (e.g. compare the drop by 
a factor of 11 at / = 0.05 Hz, to the decrease by a factor 
of 6 at / = 2 Hz, see fig. EB). Another remarkable fea¬ 
ture is that the way the linear loss modulus varies with 
frequency in samples with NPs differs from that without 
nanoparticles. In both cases, the loss modulus decreases 
when the frequency increases, except at large frequency 
(/ larger than 1 Hz), where G” eventually reaches a 
plateau or slightly increases. However, we find that the 
evolution of G” with frequency is much smoother in the 
presence of NPs (fig. SB). 

Overall, these data show that the behavior of the loss 
modulus is rather complex. Naively, one would expect 
that dissipation increases as the crystal contains more 
defects. When the heating rate increases, the crystal¬ 
lization process is faster, yielding smaller crystallites and 
more defects. Thus, one would expect larger G” val¬ 
ues as T increases. For samples without NPs, we indeed 
measure a continuous increase of G” with T (circles in 
fig. EB). Although the effect is rather weak and only 
detectable at high frequency, the measured trend is in 
agreement with the physical expectation. The results for 
samples with NPs are in sharp contrast with this picture. 
We find that dissipation is strongly influenced by the 
way the samples are prepared. Surprisingly, we observe 
a marked decrease of G” when T increases (squares in 
fig-EB), a trend opposite to that measured without NPs. 




Figure 3: (Color online). Elastic plateau (A) and loss modu¬ 
lus (B), as a function of the temperature ramps used to pre¬ 
pare the samples. Samples without nanoparticles and with 
various amounts of nanoparticles, as indicated in the caption, 
are compared. The loss modulus has been measured at 2 Hz. 
The lines are guides for the eye. 


Microscopy and scattering experiments 0113 indicate 
that when T increases the average size of the crystallite 
decreases. Hence, the total surface of grain boundaries, 
where we expect dissipation to occur predominantly, in¬ 
creases. We are then left with the puzzling observation of 
a growth of the regions where dissipative processes tend 
to concentrate, concomitant with a decrease of the loss 
modulus. 


Correlation with structural data 

These counterintuitive results can be rationalized by 
considering the structure of the samples at length scales 
comparable to the size of the nanoparticles. Figure [5] 
shows the scattering intensity profiles of several samples 
with and without NPs, comprising different mixtures of 
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Figure 4: (Color online). Strain sweeps (A) and frequency 
sweeps (B) for samples with 1 % nanoparticles, prepared using 
various temperature ramps, whose rate T is indicated in the 
legend (symbols). / = 0.5 Hz in the strain sweep experiments 
and 7 = 0.1% in the frequency sweep experiments. In (B) only 
the loss modulus is plotted. The continuous lines correspond 
to the data for a sample without nanoparticles, and with T = 
2 °C/min. 

H 2 O and D 20 ^ but all prepared with a same heating rate 
(T = 0.007 °C/min). The sample without NPs exhibits 
structural peaks due to the crystalline order of the mi¬ 
celles in the material (green line and arrows in fig. Ej). 
The same peaks are also found for samples with NPs, 
thus demonstrating that NPs do not modify the crys¬ 
talline packing (blue circles). In the presence of NPs, an 
additional peak located at g* = 0.126 nm“^ is detected. 
In experimental conditions where the scattering length 
density of the polymer is matched to that of the solvent, 
the signal of the colloidal crystal vanishes, and only the 
signal from the NPs is visible. In this case, the peak 
at g* is still present, demonstrating that it originates 
from the structure of the NPs in the colloidal crystal. 
The peak position yields an average distance between 
the NPs, 27r/g* = 50 nm. This length scale is slightly 
larger than the NP diameter (30 nm) and corresponds to 


an effective volume fraction of the NPs of the order of 
(10 — 15)%, hence a concentration very large compared 
to the average concentration of NPs in the sample (1%). 
We interpret this peak as the signature of the packing 
of the NPs in the grain boundaries. The position of this 
peak does not vary significantly with 0 or T, indicating 
a nearly constant average distance between NPs in the 
GBs [291 . By contrast, its amplitude changes. Note that 
the amplitude of the peak is only dictated by the total 
number of NPs in the GBs, as the NPs in the bulk of 
the crystalline grains are too diluted to contribute to the 
scattering intensity in the range of wave vectors inves¬ 
tigated here. We quantify the relative amplitude of the 
peak by ^peak = {I* - Imin)/Iinin, where I* is the scat- 
tered intensity measured at the peak position, g*, and 
/min is the minimum scattered intensity measured at g 
lower than g*, and plot the evolution of Apeak with T 
and in figs. [6B and D. For samples comprising a fixed 
amount of nanoparticles, we find that Apeak decreases 
as T increases (fig. [6^). Our data thus indicate that for 
these samples the total amount of NPs in the GBs de¬ 
creases when T increases. However, as T increases the 
average crystallite size decreases, and thus the total sur¬ 
face area of GBs increases. The fact that the cumulated 
amount of NPs in the GBs decreases, in spite of an in¬ 
crease of the total surface area of the grain boundaries, is 
due to a change in the efficiency of the segregation pro¬ 
cess of the NPs during crystallization. Indeed, as already 
briefly mentioned in the Materials and Methods section 
and as described in detail elsewhere [23|, not all the NPs 
incorporated in the polycrystal are located in the grain 
boundaries: the partitioning coefficient V depends on the 
sample microstructure. For fluorescent polystyrene NPs 
of size comparable to the silica NPs used here, we have 
previously found that V diminishes as the average size of 
the crystallites decreases. Although the NPs investigated 
here are of different chemical nature, it is likely that a 
similar behavior also holds for the silica NPs. We indeed 
found that the evolution of the grain size with (j) and T 
are similar for both kinds of particles, although the grains 
are slightly larger for silica NPs than for polystyrene 
NPs. Additionally, a similar dependence on crystalliza¬ 
tion rates of the partitioning coefficient is commonly ob¬ 
served in atomic and molecular systems [s^. The data 
shown in fig[6B suggest that, in the linear regime, dissi¬ 
pation is mostly due to the NPs that are segregated in 
the grain boundaries. Within this scenario, the decrease 
of G” with T would be the consequence of reduced dis¬ 
sipation stemming from a smaller number of NPs in the 
grain boundaries as T increases, due to the reduction of 
the partitioning coefficient. 

To further test these ideas, we perform rheology and 
SANS measurements on a series of samples with various 
amounts of NPs, keeping fixed the heating rate. Data 
obtained by SANS (fig. [6p) show a non-monotonic be¬ 
havior of Apeak with (j). This can be interpreted as the 










6 



q (nm^) 



0.001 0.01 0.1 


r(°C/min) 


0.0 0.5 1.0 1.5 2.0 


Figure 5: (Color online). Small-angle neutron scattering in¬ 
tensity profiles for samples prepared with a temperature ramp 
of T = 0.007 °C/min. The intensity is normalized by the 
scattering contrast Ap between the solvent and the scatter¬ 
ing objects (polymer, blue circles, and NPs, red squares). The 
green continuous line is a sample without nanoparticles (NPs) 
and the symbols correspond to samples comprising 1 % NPs 
under different scattering contrast conditions. Red squares: 
polymer-matched sample, only the micelles contribute to I{q). 
Blue circles: silica-matched sample. The dotted lines indicate 
the position of the structure peak of the NPs, at g*, due to 
their confinement in the grain-boundaries, and the intensity 
/* and /min used to measure the peak amplitude. The ar¬ 
rows indicate the peaks due to the scattering of the micellar 
crystal. 


result of the competition between two opposite effects: 
on the one hand, the higher 0, the larger the number 
of NPs that may be segregated in the GBs. On the 
other hand, the higher (j) the lower the partition coef¬ 
ficient and hence the lower the fraction of NPs that ac¬ 
tually accumulate in the GBs rather than in the bulk of 
the crystallites. At small 0, the first effect dominates: 
the size of the grains and the partition coefficient V do 
not vary significantly with 0, implying that the cumu¬ 
lated number of NPs in the GBs is directly proportional 
to the concentration of NPs. Gonsequently, Apeak ^ 0, 
consistently with the experimental data (dotted line in 
fig- EP)- As (j) increases further, the average size of the 
crystallite decreases and the partitioning becomes less 
efficient, leading to a decrease of the relative amount of 
NPs in the grain-boundaries. This explains why Apeak 
stops growing and eventually slightly decreases at higher 
(j). Remarkably, a similar non-monotonic behavior with 
(j) is observed also for the loss modulus (fig. [6p), lead¬ 
ing support to the hypothesis that G” is dominated by 
the dissipation associated with the nanoparticles in the 
grain-boundaries. 


Figure 6: (Color online). Comparison between the lin¬ 
ear visco-elasticity measurement and the structural measure¬ 
ments. In (A,C), the loss modulus taken at 2 Hz is plotted as 
a function of the heating ramp used to prepared the samples, 
T (A), and as a function of the concentration of nanoparticles, 
(j) (C). In (B,D), the relative amplitude of the NPs structure 
peak is plotted as a function of T (B), and as a function of 
0 (D). In (A,B) 0=1%. In (C), T = 0.001 °C/min, and in 
(D), T = 0.02°C/min. In (D), the dotted line is a guide for 
the eye. 

Discussion 

We have found that the linear storage modulus is 
nearly frequency independent and constant irrespective 
of the presence or not of nanoparticles and of the heat¬ 
ing rate used to prepare the samples, within experimental 
uncertainties. An average over all measurements yields 
Go = (11970 ± 750) Pa. Consistently with previous re¬ 
sults obtained in comparable conditions 0 , we do not 
measure any weakening of the elasticity due to the pres¬ 
ence of NPs. The elasticity is independent of the presence 
of NPs and of the heating rate, and thus on the crystal 
microstructure. Hence, the elasticity is governed solely 
by the packing of the micelles in the crystalline lattice. 
Using simple scaling arguments. Go is expected to scale 
as ksT/a^^ where a = 11 nm is the radius of the micelles. 
We find Goa^{ksT = 3.9, a value in good agreement with 
experiments on colloidal hard sphere crystals, for which, 
e.g., Goa^lksT = 5.8 at a colloid volume fraction of 
0.55 [13. 

The behavior of the loss modulus is much richer. On 
the one hand, the physical origin of the dissipation with 
and without nanoparticles is different, as reflected by the 
very different dependence of G” with frequency in the 
two cases. For samples without NPs, the large values 
of the loss modulus at small frequency indicate dissipa¬ 
tion processes with long characteristic relaxation times, 
presumably due to point and line defects in the bulk of 
the crystallites, since G" is insensitive to the amount of 










GBs at low frequency. Samples with NPs and prepared 
at a high heating rate certainly contain a large number 
of bulk defects, since the NP impurities tend to be dis¬ 
persed throughout the sample, rather than massively ac¬ 
cumulate in the GBs. One might then expect that G" 
would be as high or even higher than in the NPs-free 
samples. Surprisingly, in the presence of nanoparticles 
the opposite trend is observed: G" has lower values and 
is almost frequency-independent. We propose that this 
stems from a stabilizing effects of the NPs, leading to 
a slowing down of the defect dynamics, as observed in 
colloidal crystals [13, [HI and in other materials jHI-[4l| . 

For samples with NPs but prepared at lower heating 
rates, the concordance between the SANS and the visco¬ 
elasticity measurements suggests that the loss modulus is 
directly related to the total amount of NPs in the GBs, 
and thus that dissipation processes mainly involve the 
NPs confined in the grain boundaries. These processes 
are likely to be influenced by any interaction between 
the polymer and the NPs, such that the surface chem¬ 
istry of the NPs should play a role. For the silica NPs 
used here, it is known that the copolymer FI08 slightly 
adsorbs to the particle surface [29|, 1^. Note, in addi¬ 
tion, that the fact that the structure peak originating 
from the NPs accumulated in the GBs is also visible 
for a SANS profile measured in silica-matched conditions 
(red squares in fig. [5]), is also an indication of interac¬ 
tions between the F108 polymer and the NPs. We test 
the effect of the NPs-F108 interactions by performing 
a few additional visco-elasticity measurements using a 
different kind of nanoparticles, carboxylated polystyrene 
latexes, of diameter 36 nm, comparable to that of the 
silica NPs. The effect of the polystyrene NPs on the 
sample structure, partial accumulation in the GBs and 
influence on the microstructure, is similar to that of the 
silica NPs [23|. Moreover, experiments on dilute suspen¬ 
sions of polymer and polystyrene NPs indicate an adsorp¬ 
tion of the central block of the polymer to the surface 
of the NPs [i^ with an amount of adsorbed polymer 
per unit area of NPs coinparable to the one evaluated 
for silica NPs (Img/m^) [HI- Despite these similarities, 
the impact of polystyrene NPs on the sample dissipation 
is strikingly different to that on silica NPs. This is il¬ 
lustrated in fig. 13 where the frequency-dependent loss 
modulus is plotted for two heating rates T = l°C/min 
and T = 0.01°G/min, for the two kinds of NPs. The 
significant increase of the loss modulus at low / observed 
for silica NPs-laden samples prepared at a small heating 
rate is not found for samples comprising polystyrene NPs. 
By contrast, we find that, for a sample with polystyrene 
NPs, G” is comparable for samples prepared with a fast 
heating rate and with a slow heating rate. Hence at 
high rates, the behavior of samples comprising silica or 
polystyrene NPs are comparable, whereas they markedly 
differ at a low rate. This is consistent with the gen¬ 
eral picture that at low rate dissipation occurs mainly 



Figure 7: (Color online). Frequency dependence of the 
loss modulus for samples with nanoparticles, prepared us¬ 
ing two different temperature ramps, l°C/min (circles) and 
0.001 °C/min (squares). Two kinds of particles are used: sil¬ 
ica NPs (red, dark gray, symbols) and polystyrene (PS) (yel¬ 
low, light gray, symbols). The volume fraction of NPs is 
(j) = 0.5 %. The strain amplitude is 0.1%. For the sake of 
comparison, the black continuous line shows the data for a 
sample without nanoparticles prepared with T — 1 °C/min. 

in the GBs (where NPs-polymer interactions are impor¬ 
tant), while at high rate NPs are more evenly dispersed 
in the sample and act as point defects that stabilize the 
sample. These preliminary, intriguing results deserve a 
deeper investigation that is out of the scope of the present 
work. They nevertheless illustrate nicely the crucial role 
of the interaction between the polymeric micelles and the 
NPs. 

CONCLUSION 

We have investigated how impurity-like nanoparti¬ 
cles added in a colloidal crystal impact its linear visco¬ 
elasticity, because of their partial segregation in the 
grain-boundaries that separate crystalline grains with 
different orientations. Two independent experimental 
quantities, the volume fraction of NPs and the heating 
rate that controls the crystallization, have been varied, 
both allowing the microstructure of the polycrystal to 
be tuned. We have found that the elasticity is inde¬ 
pendent of the presence of nanoparticles and the aver¬ 
age size of the crystallites. In sharp contrast, the dis¬ 
sipation strongly varies with the presence of NPs and 
with the sample microstructure. By combining rheology 
measurements and structural measurements, we conclude 
that the viscous dissipation is governed by the cumulated 
amount of nanoparticles located in the grain boundaries. 
We are not aware of previous comparable results for col¬ 
loidal crystals. A comparison can nevertheless be drawn 
with other soft materials. How the properties of the in- 
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terface impact the loss modulus has for instance been 
investigated in foams Q. Interestingly, in particles- 
laden foams the loss modulus G” has been observed to 
increase with the amount of particles [45|, indicating, as 
in our case, that dissipation processes are correlated to 
the presence of nanoparticles in the interface. However, 
in particles-laden foams the storage modulus increases 
as well with the amount of particles, because the elas¬ 
ticity of foams is also directly related to the interfacial 
properties. By contrast, in our poly crystals elasticity 
and dissipation are governed by different structural fea¬ 
tures, the packing of the micelles for the former, and 
the amount of NPs in the grain boundaries for the lat¬ 
ter. In this respect, colloidal polycrystals doped with 
impurities share more analogies with their atomic and 
molecular counterparts. The fact that impurity segre¬ 
gation in the grain-boundaries modifies the mechanical 
properties of the polycrystal is indeed reminiscent of the 
embrittlement of metal by solute segregation in the grain 
boundaries that is eventually responsible for intergranu¬ 
lar fracture leading to material failure Hi-lii. 


ACKNOWLEDGEMENT 

This work has been supported by ANR under Contract 
No. ANR-09-BLAN-0198 (COMET), and by Laboratoire 
Leon Brillouin (LLB). We thank F. Cousin (LLB) for help 
with the small-angle neutron scattering measurements. 
Discussion with N. Ghofraniha, F. Caton and D. Roux 
are gratefully acknowledged. 


* Electronic address: laurence.ramos@umontpellier.fr 

[1] J. M. Howe, Interfaces in Materials, (Wiley, New-York, 
1997). 

[2] L. Lu, M. L. Sui, and K. Lu, Science, 287, 1463, 200. 

[3] J. Schi0tz, F. D. Di Tolla, and K. W. Jacobsen, Nature, 
391, 561, 1998. 

[4] J. Schi0tz, and K. W. Jacobsen, Science, 301, 1357, 2003. 

[5] V. Yamakov, D. Wolf, S. R. Phillpot, A. K. Mukherjee, 
and H. Gleiter, Nature Materials, 3, 43, 2004. 

[6] D. Pan, S. Kuwano, T. Fujita, and M. W. Chen, Nanolet¬ 
ters, 7, 2108, 2007. 

[7] H. Li, H. Choo, Y. Ren, T. A. Saleh, U. Lienert, P. K. 
Liaw, and F. Ebrahimi, Phys. Rev. Lett., 101, 015502, 
2008. 

[8] Y. M. Wang, R. T. Ott, A. V. Hamza, M. F. Besser, J. 
Aimer, and M. J. Kramer, Phys. Rev. Lett., 105, 215502, 
2010. 

[9] T. Chookajorn, H. A. Murdoch, and C. A. Shuh, Science, 
337, 951, 2012. 

[10] Z. Shan, E. A. Stach, J. M. K. Wiezorek, J. A. Knapp, D. 
M. Follstaedt, and S. X. Mao, Science, 305, 654, 2004. 

[11] J. P. Buban, K. Matsunaga, J. Chen, N. Shibata, W. Y. 
Ching, T. Yamamoto, and Y. Ikuhara, Science, 311, 212, 
2006. 


[12] S. Cheng, A. D. Stoica, X.-L. Wang, Y. Ren, J. Aimer, 
J. A. Horton, C. T. Liu, B. Clausen, D. W. Brown, P. K. 
Liaw, and L. Zuo, Phys. Rev. Lett., 103, 035502, 2009. 

[13] M. Zapotocky, L. Ramos, P. Poulin, T. C. Lubensky,and 

D. A. Weitz, Science, 283, 209, 1999. 

[14] T. A. Wood, J. S. Lintuvuori, A. B. Schofield, D. Maren- 
duzzo, and W. C. K. Poon, Science, 334, 79, 2011. 

[15] E. Eiser, F. Molino, and C. Porte, Eur. Phys. J. E 2 , 39 
( 2000 ). 

[16] E. Eiser, F. Molino, C. Porte, and X. Pithon, Rheol. 
Acta, 39, 201, 2000. 

[17] D. C. Pozzo, K. R. Hollabaugh, and L. M. Walker, J. 
Rheol. 49, 759 (2005). 

[18] M. A. Torija, S.-H. Choi, T. P. Lodge, and F. S. Bates, 

J. Phys. Chem. B, 115, 5840, 2011. 

[19] T. Palberg, W. Monch, J. Schwarz, and P. Leiderer, J. 
Chem. Phys., 102 , 5082, 1995. 

[20] A. Engelbrecht, and H. J. Schope, Crystal Growth & 
Design, 10, 2258, 2010. 

[21] S. Gokhale, K. H. Nagamanasa, R. Ganapathy, and A. 

K. Sood, Soft Matter, 9, 6634, 2013. 

[22] Y. Peng, Z. Wang, A. M. Alsayed, A. G. Yodh, and Y. 
Han, Phys. Rev. Lett., 104, 205703, 2010. 

[23] N. Ghofraniha, E. Tamborini, J. Oberdisse, L. Cipelletti, 
and L. Ramos, Soft Matter 8, 6214 (2012). 

[24] S. Gokhale, K. H. Nagamanasa, V. Santosh, A. K. Sood, 
and R. Ganapathy, Proc. National Acad. Sci. USA, 109, 
20314, 2012. 

[25] P. Lejcek, Grain Boundary Segregation in Metals 
(Springer, 2010). 

[26] W. Losert, B. Q. Shi, and H. Z. Cummins, Proc. National 
Acad. Sci. USA, 95, 431, 1998. 

[27] V. W. A. de Villeneuve, R. P. A. Dullens, D. G. A. Aartz, 

E. Groeneveld, J. H. Scherff, W. K. Kegel, and H. N. K. 
Lekkerkerker, Science, 309, 1231, 2005. 

[28] K. Yoshizawa, T. Okuzono, T. Koga, T. Taniji, and J. 
Yamanaka, Langmuir, 27, 13420, 2011. 

[29] E. Tamborini, N. Ghofraniha, J. Oberdisse, L. Cipelletti, 
and L. Ramos, Langmuir 28, 8562 (2012). 

[30] A. Louhichi, E. Tamborini, N. Ghofraniha, F. Caton, 

D. Roux, J. Oberdisse, L. Cipelletti, and L. Ramos, Phys. 
Rev. E 87, 032306 (2013). 

[31] E. Tamborini, L. Cipelletti, and L. Ramos, Phys. Rev. 
Lett., 113, 078301, 2014. 

[32] P. Alexandridis, T. Nivaggioli, and T. A. Hatton, Lang¬ 
muir, 11 , 1468, 1995. 

[33] R. G. larson. The Structure and Rheology of Complex 
Fluids, (Oxford University Press, New York, 1999). 

[34] P. Lindner, X-ray and Light Scattering, North Holland, 
Elsevier, 2002. 

[35] H. M. Wyss, K. Miyazaki, J. Mattsson, Z. Hu, D. R. Re- 
ichman, and D. A. Weitz, Phys. Rev. Lett., 98, 238303, 
2007. 

[36] M. J. Aziz, J. Y. Tsao, M. O. Thompson, P. S. Peercy, 
and C. W. White, Phys. Rev. Lett., 56, 2489, 1986. 

[37] S.-E. Phan, M. Li, W. B. Russel, J. Zhu, P. M. Chaikin, 
and C. T. Lant, Phys. Rev. E, 60, 1988, 1999. 

[38] V. W. A. de Villeneuve, L. Derendorp, D. Verboekend, 

E. C. M. Vermolen, W. K. Kegel, H. N. K. Lekkerkerker, 
and R. P. A. Dullens, Soft Matter, 5, 2448, 2009. 

[39] W. A. Soer, J. Th. M. De Hosson, A. M. Minor, J. W. 
Morris Jr., and E. A. Stach, Acta Materialia, 52, 5783, 
2004. 

[40] Y. Li, D. Raabe, M. Herbig, P.-P. Choi, S. Goto, A. 




9 


Kostka, H. Yarita, C. Borchers, and R. Kirchheim, Phys. 
Rev. Lett., 113, 106104, 2014. 

[41] P. M. Hazzledine, and R. D. J. Oldershaw, Phil. Mag. A, 
61, 579, 1990. 

[42] J. A. Shar, T. M. Obey, and T. Cosgrove, Colloids Surf. 
A, 150, 15, 1999. 

[43] J.-T. Li, K. D. Caldwell, and N. Rapoport, Langmuir, 
10, 4475, 1994. 

[44] S. Marze, R. M. Guillermic, and A. Saint-Jaimes, Soft 
Matter, 5, 1937, 2009. 


[45] S. Cohen-Addad, M. Krzan, R. Hohler, and B. Herzhaft, 
Phys. Rev. Lett., 99, 168001, 2007. 

[46] G. Duscher, M. F. Chisholm, U. Alber, and M. Riihle, 
Nature Materials, 3, 621, 2004. 

[47] R. Schweinfest, A. T. Paxton, and Michael W. Finnis, 
Nature, 432, 1008, 2004. 

[48] H.-P. Chen, R. K. Kalia, E. Kaxiras, G. Lu, A. Nakano, 
K. 1. Nomura, A. C. T. van Duin, P. Vashishta, and Z. 
Yuan, Phys. Rev. Lett., 104, 155502, 2010. 


